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Abstract 

We investigate the prospects for the discovery of a neutral Higgs boson produced in association 
with a b quark, followed by the Higgs decay into a pair of bottom quarks, pp — > b(j)° — > bbb + 
X, at the CERN Large Hadron Collider (LHC) within the framework of unified supersymmetric 
models. The Higgs boson can be a heavy scalar H° or a pseudoscalar ^4°. Furthermore, this 
direct discovery channel is compared with the indirect Higgs searches in the rare decay B s — > 
at hadron colliders. Promising results are found for the minimal supergravity (mSUGRA) 
model, the anomaly mediated supersymmetry breaking (AMSB) model, and the gauge mediated 
supersymmetry breaking (GMSB) model. We find that the indirect search for B(B S — > u + fi~) > 
5 x 10~ 9 is complementary to the direct search for b(p° — > bbb with yfs = 14 TeV and an integrated 
luminosity (L) of 300 fb _1 . In the AMSB and GMSB models, b(f>° -> bbb with L = 300 ftr 1 covers 
a larger area in the parameter space than B(B S — > u + /j,~) > 5 x 10~ 9 . In addition, we present 
constraints from b — > s'j and muon anomalous dipole moment (Aa^) on the parameter space. 
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I. INTRODUCTION 



Although the Standard Model (SM) remains an incredibly successful description of most 
of the phenomena studied by high energy physics, it is expected to be replaced by a more 
complete model in order to address a number of theoretical problems. There are also several 
experimental signatures for which the Standard Model does not seem to account. Perhaps 
the most pressing theoretical question is the stabilization of the hierarchy between the weak 
and Planck scales. Supersymmetry (SUSY) provides the most popular solution to this 
problem and, due to the large number of exotic particles introduced, also suggests solutions 
to a variety of other Standard Model anomalies. These include the origin of dark matter 
and deviations in flavor and CP physics. On the other hand, SUSY must be a broken 
symmetry at low scales and in the absence of a detailed theory of spontaneous breaking we 
are left with a large number of free parameters in the form of soft supersymmetry breaking 
terms. To ameliorate this profusion of parameters, it is useful to adopt one of several unified 
frameworks which reduce the number of new variables. 

For phenomenological reasons, it is assumed that supersymmetry breaking is driven by 
the dynamics of a hidden sector separate from the Standard Model, and this breaking is com- 
municated to the Standard Model particles and superpartners by a messenger sector which 
couples to both. Gravity will act as a natural candidate for this messenger, leading to super- 
gravity theories (SUGRA). However, it is possible that the gravity-induced SUSY breaking 
terms are not dominant and that the leading effects are generated either by messenger fields 
which carry Standard Model gauge numbers (gauge-mediated supersymmetry breaking, or 
GMSB), or by terms arising from the superconformal anomaly (anomaly- mediated super- 
symmetry breaking, or AMSB). 

It is also a striking fact that adding the minimal particle content required to incorporate 
the SM fields in a supersymmetric theory changes the running of the gauge couplings such 
that they unify at a high scale (approximately 10 16 GeV). This is a much better unification 
than is seen for the SM fields alone and is consistent with the Grand Unified Theory (GUT) 
hypothesis that the three Standard Model gauge groups are remnants of a larger symmetry 
which is restored at high scales. 

With the assumptions of grand unification and the messenger sector for SUSY breaking, 
the theory can be characterized by measured Standard Model inputs and a small number 
of additional unification scale parameters. This is enough to determine all the masses and 
couplings of all the superpartners, allowing us to predict the potential for direct discovery 
and to calculate indirect effects, which appear as corrections to precision measurements 
relative to their Standard Model values. 

The LHC is rapidly accumulating data, currently with y/s = 7 TeV, and already setting 
new limits on observables sensitive to new physics. It is planned to run at this energy 
through 2012, followed by a long shutdown for refurbishing and then a new run in 2014 
aiming at the design energy of \fs = 14 TeV. With high energy and luminosity, the LHC 
should be able to discover both the Higgs boson and supersymmetry if they exist in the 
100 — 1000 GeV range as anticipated by theorists. 

One of the basic predictions of supersymmetry is the existence of (at least) two Higgs 
doublets, in contrast with the Standard Model. This leads to five physical states at low 
energy: two neutral scalar s h° (lighter) and H° (heavier), a neutral pseudoscalar A , and a 
pair of charged scalars H . Thus a key signature of supersymmetry would be the detection 
of this extended Higgs sector. In this paper, we focus on the prospects for detecting the 
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heavy neutral Higgs scalar and the Higgs pseudoscalar. We investigate the decay channel 
0° — > bb,(j) = H°,A°, which becomes important for large values of the ratio of vacuum 
expectation values (VEVs) (tan/3 = v^/vi). 

We have considered this channel in a previous paper in the context of the minimal su- 
persymmetric standard model (MSSM) pj. In this study, we choose instead to work in the 
three unifying frameworks mentioned above. We plot the estimated detection contour in the 
space of the unification scale parameters for mSUGRA, mGMSB, and mAMSB respectively. 
These models also make predictions for rare decay rates and precision parameters which are 
currently measured or constrained. In particular, we compare the LHC discovery potential of 
the 36 channel with B s — > fi + at hadron colliders, which is especially enhanced by a large 
value of tan /3. Furthermore, we consider constraints from b — > s'j and and the anomalous 
magnetic moment of the muon, (g — 2) M . 

In our mSUGRA and mAMSB plots, the lightest neutralino \ is always required to be the 
lightest SUSY particle, or LSP. 1 If i?-parity is conserved, then the x will De absolutely stable 
and may comprise all or some of the cold dark matter (CDM) in the universe. Within the 
confines of the MSSM, a neutralino as CDM requires rather high fine-tuning [2] and also very 
low reheat temperature Tr < 10 5 GeV in order to avoid bounds on late-decaying gravitinos 
from BBNjJjj]. If we invoke late-decaying scalar fields, such as TeV-scale moduli, into our 
theory, then these can either increase or decrease the standard neutralino abundance^. 
Alternatively, invoking the Peccei-Quinn solution to the strong CP problem then requires 
either mixed axion/axinojBf or mixed axion/neutralino CDMpJ, for which again the relic 
density predictions can be quite different. Here, we will simply refrain from implementing 
any such constraints. 



II. HIGGS PRODUCTION IN ASSOCIATION WITH B QUARKS 

The Minimal Supersymmetric extension of the Standard Model (MSSM) requires two 
Higgs doublets with opposite hypercharge to account for fermion masses since the superpo- 
tential must be analytic: 



'J \ _ I ^2 



(1) 



Each doublet acquires a VEV in its neutral component, V\ and The model is char- 
acterized by the ratio of these VEVs tan/3 = v 2 /vi. After symmetry breaking three of the 
degrees of freedom between the doublets are eaten by the massive weak vector bosons and 
the remaining five become the massive physical eigenstates h°, H°, A and if ± which are a 
mixture of the original doublet fields. By convention, h° is lighter than H°. We will use 0° 
to refer generically to any of the neutral Higgs bosons. 

At the LHC, 0° is predominantly produced by gluon fusion (gg — > 0°) for tan/3 < 5. 
However, in a 2HDM the coupling of down-type quarks to the Higgs is proportional to the 
quark mass divided by cos /3, which can be quite large compared to the Standard Model value. 
Hence for tan/3 > 7 the leading production of </>° comes from bb fusion jZ-U|. Similarly, the 
branching fraction for 0° — > bb can become nearly 90%. 



1 In mGMSB, the gravitino G is taken as LSP. 
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In the case of high tan f3, the simplest channel to look for neutral Higgs would seem to be 
bb — > — > bb. However, this signal is swamped by large QCD dijet backgrounds. The next 
step is to take advantage of b-tagging capabilities by looking for 0° in association with one 
or more additional b-jets. Here we have the option of requiring four b-tagged jets with high 
Pti PP <j>°bb —> bbbb + X, which will greatly reduce the background compared to the dijet 



case but which suffers from a small signal in turn [12|-|l4j. In this paper, we will pursue the 
intermediate case with 3 high-p^ tagged b- jets in the final state. It has been argued that 
this channel is more promising at the LHC 
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We treat this channel in a 5-flavor quark Parton Distribution Function (PDF) scheme. In 
a 4-flavor scheme, the leading order process of interest would be gg — > <f)°bb, where the Higgs 
subsequently decays back to b-jets. One would then need to sum over all configurations 
with 3 high-pr quarks in the final state. This requires careful treatment of higher-order 
corrections since the 4th quark, at low p?, gives rise to large leading- log corrections. In 
the 5-flavor scheme the b-quark is treated as an initial parton and the large corrections are 
absorbed by the b-quark PDF. For 5 flavors the leading order process is bg — > bcfp — > bbb. 
(As well as the conjugate process with 60°. Henceforth we will use b to refer to both b 
and anti-6 quarks.) The process bg — > b<p° has two Feynman graphs and by choosing the 
factorization and renormalization scale /ip = = Mfl-/4, NLO corrections can be made 



relatively small [10|, [15 



The decay of the Higgs boson 0° — > bb depends at first order on the size of the Yukawa 
coupling to bottom quarks, which scales with tan/3. Here, we take into account several 
higher order effects since we wish the effective size of this vertex to correspond with the 



Higgs decay width to bb. The decay width at NLO can be written |l_6 



k0 ~* hl) = 3 ~^r mlm[1 + Aqcd + At]{gt)2 > (2) 



where 



A QCD = 5.67^^ + (35.94 - 1.36iV F )(^^) 2 (3) 

7T 7T 



A^cot^[3.83-ln(M ) + l ln( ^ ))2]( ^) )2 (4) 
A--c t«c t /3[ 1.57^1n4 + il n 4 ) 1(^ 



(5) 



Af = -tanacot/3[1.57- -ln(^) + - l n (^)) 2 ](^— ^) 2 (6) 



H u lv± t 
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— sin a A 



cos /3(1 + A b ) tan a tan /3 



1 " I f— o) (7) 



a H _ cos a A b 
yb cos/3(l + A 6 ) 1 cotatan/T 1 ; 



rf = i^-d - r% 0) 

1 + A& tan z p 

The quantity A;, is an effective shift in the b-quark mass induced by SUSY QCD 
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corrections 17H20j. It is computed from the sum of two terms, = + A£. 

2a 



A b = ~^ m 9^ taI1 P I ( rn b 1 > m b2> rn 9) ( 10 ) 



= — A t fi tan (31 (m h ,m i2 ,fi) (11) 

6 ' C) = ' (a 2 -^)(^- c2 ) (c2 - Q2 ) ( QV ln ^ + &V L 4 + C2fl2 111 ^ (12) 

In the equations above a t = At/(47r), where At = y2m t /v2 is the top Yukawa coupling, and 
A t is the trilinear stop coupling. A b can be a significant correction to the decay rate. One 
can see that Ajj decreases the rate for a positive sign of \i and increases the rate for negative 
/i. The A^ term may be comparable in size or quite small, with a variable sign depending 
on the trilinear coupling. In addition to the Higgs decay vertex, we include a At correction 
factor at the production vertex. 

In order to evaluate the expressions above, as well as the Higgs mass itself, we use the 
program ISAJET, which implements one of several available SUSY-breaking schemes in 
terms of GUT scale input parameters, then computes the low scale SUSY masses based on 
running of the parameters according to the renormalization group equations j2lj . We use the 
full width of the Higgs as computed by ISAJET, modified by the corrected bb decay width 
as indicated above. In general, the full width of the Higgs is small enough in comparison 
to the Higgs mass to use the Narrow Width Approximation (NWA)[lj. For large masses 
and high tan/3 the NWA becomes slightly worse. In practice we use a full Breit-Wigner 
calculation for bg — > 60° — » bbb with the coupling of the two outgoing 6s from Higgs decay 
set to an effective value which reproduces the bb decay width. 



We simulate the signal using MadGraph4 [22|-|24j to generate our matrix elements (bg — >■ 
bbb) and then evaluate the total cross-section with a Monte Carlo program. We smear the 
momenta of the outgoing b-quarks with a Gaussian distribution parametrized by 

AE 0.60 , , 

— = 7f ®0.03, (13) 



based on ATLAS estimates of detector effects [25f| . We then impose the following cuts on 
our signal: 

• pt > 70 for all three jets, 

• \tj\ < 2.5 for all three jets, 

• ARij > 0.7 for each pair of jets (i,j), 

• Missing E T < 40 GeV, 

• | My — M^| < 0.15M<£ for at least one pair of jets. 

In the list above r/ is the pseudorapidity and AR = a/ Arj 2 + A</> 2 is a measure of jet 
separation. We choose the p?, rj, and AR cuts to provide for good reconstruction of three 
jets in the b-tagging region. The px cut is also chosen in accord with the CMS Level 1 



Trigger for 3 jets |26J. Mj,- is the invariant mass constructed from two of the three jets; we 



require that it be within a 15% window of the true pseudoscalar mass Ma for at least one 
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pair of jets. We found in a previous work that allowing only the two highest px jets for this 
requirement produces a slight improvement with respect to background at high masses 

The backgrounds to our signal are dominated by pure QCD processes which produce three 
jets. The irreducible background is bg — > bbb. We also include the reducible backgrounds: 
pp — > ebb + X, pp — > gbb + X, and pp — > qbb + X where q = u,d,s along with pp — >■ tt — >■ 
bbjjlv + X and pp — > tt — >■ bbjjjj + X. These involve one mis-tagged non-b jet and can 
potentially be reduced with improved b-tagging. For both the signal and the background 
we assume an effective b-tag rate of €b = 0.6 with a mis-tag rate e c = 0.14 for c-quark jets 
and €j = 0.01 for light jets (g,u,d,s). With these efficiencies, the pure QCD backgrounds 
(bbb, ebb, gbb, qbb) are all of comparable size, with bbb or gbb being the largest, while the 
backgrounds with intermediate t-quarks are negligible in comparison. We include them 
nonetheless for completeness. The double mis-tagged background ccb should be roughly ^ 
the size of the bbb background after tagging and thus only a 1 — 2% correction to the total 
background. 

For simulation, backgrounds are generated using MadGraph amplitudes with the renor- 
malization and factorization scales set to p^(l)/2, half the transverse momentum of the 
leading jet, for the pure QCD backgrounds, and y/s for the tt backgrounds. Cuts are the 
same as used for the signal, except that we also apply a veto to tt events with 4 jets having 
Pt > 15 GeV. We assume a K factor of two for each background while keeping the K factor 



for the signal at one [12|, |27 



Once we have both signal and background we are able to compute a statistical significance, 
for which we use 

Nq 

Nss = iVs (14) 
VN S + N B 1 } 

where N$ and Nb are the number of expected signal and background events respectively. 
We set Nss = 5 as the discovery limit. In practice, we use the following process to find the 
discovery contour: First we generate a set of background estimates over a range of masses 
Ma- The background cross-sections only depend on the mass through the location and 
size of the invariant dijet mass cut as given above. Next we choose a point in the GUT- 
scale parameter space for one of the SUSY-breaking models and use Isajet to calculate 
the relevant weak scale parameters, namely, the Higgs mass and decay width. We modify 
the (fP — > bb width with the corrections listed above and feed these parameters into our 
MadGraph/Monte Carlo program which calculates our signal for that point in parameter 
space. The background at that point is determined by cubic spline interpolation from our 
array of Ma dependent backgrounds and a significance is calculated. Then by scanning over 
one of the GUT-scale parameters with the others fixed we can locate the discovery contour 
where Nss — 5. 



III. EXPERIMENTAL CONSTRAINTS 

In addition to the estimated range of sensitivity in the 3b channel, it is interesting to 
consider other effects of the high tan /3 scenarios and model assumptions. There are a 
number of experimentally measured or constrained quantities which are quite sensitive to 
deviations from the Standard Model. Especially important at high tan /3 is the rare decay 
B s — > fi + fi~. In the Standard Model this decay is expected to have the low branching 
fraction BF(B S — > = (3.6 ±0.37) x 10~ 9 [28(. However, diagrams involving H and A 
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in a 2HDM are proportional to (tan/3) 3 , meaning that their contributions to the decay rate 
scale as (tan/3) 6 . This can change the predicted decay rate by orders of magnitude and the 
current results from several experiments put us in a very exciting time frame. 

At the Tevatron, DO and CDF set limits on the observed branching fraction of BF(B S —> 
/i + /i _ ) < 5.1 x 1CT 8 and BF(B S — > /i + /i") < 4. x 1CT 8 . Meanwhile, early results from 



the LHC have lowered this limit to 1.9 x 1(T 8 (CMS) and 1.5 x l(T 8 (LHCb)|29|, |30j, using 



approximately 300pb _1 of data. A combined analysis using CMS and LHCb data puts the 
limit at 1.1 x 10 _8 [3l|. Thus we are already constrained to limit the effects of heavy Higgs 
bosons with large tan f3, while still allowing for an enhanced branching fraction up to 3 times 
the Standard Model rate. Interestingly, CDF has reported a weak excess corresponding to 
a signal at (1.8±l#) x 1(T 8 with 98% probability of exceeding the SM rate [321] . LHCb and 
CMS have not seen evidence of this signal but have also not strongly ruled it out. LHCb in 
particular should cover the remaining space down to the Standard Model limit in the near 
future. It was anticipated to reach the SM limit with ~ 2 fb _1 of data at 14 TeV running 
(33| . Thus we expect it to show evidence of an excess signal or to put stringent limits on 
the allowed branching fraction in the next year or two. 

Supersymmetric models in general are constrained by electroweak precision data and 
searches for new flavor and CP violation. Two constraints of particular interest to us are 
the measured values of — 2 and b — » 57. Both quantities are sensitive to tan (3 (though 
not so strongly as B s — > and to the effects of the SM superpartners. 

The anomalous magnetic moment of the muon, — 2, is one of the more precisely 
calculated and measured quantities of quantum field theory. The experimental value is 
found to be = {g-2)/2 = (116592089.0 ± 6.3) x 10~ 10 [34]. A recent calculation puts the 
theoretical value for the Standard Model at a M = (11659182.8 ± 4.9) x 10~ 10 , leading to a 
discrepancy Aa M = (26.1 ±8.1) x 10~ 10 , i.e. a 3.3<r excess in experiment [35|. Other groups 
find similar results. It is tempting to attribute this excess to new physics and supersymmetric 
models can easily account for it if they have the correct set of masses and parameters. 

At high tan the dominant new contribution to a M is proportional to tan /3 with the same 
sign as \xM 2 0. (Unless \M 2 \,M% < \Mx\,M L , where Mi, Mg are the left- and right- 
handed slepton soft SUSY-breaking masses, respectively.) Thus, for the minimal models we 
consider below, one requires a positive sign for \x if MSSM contributions are to account for 
the observed excess of a M . 

A third sensitive probe is the flavor changing decay b — > sj, observed through B — > X^. 
Experimentally, this is measured at (3.55 ± 0.26) x 10~ 4 |34J|. Theoretical predictions in 
the Standard Model put the value at (3.15 ±0.23) x 10~ 4 [37|. In supersymmetric theories, 
loops involving the charged Higgs boson, as well as those involving charginos and squarks, 
can make large contributions to b — > 57. 

There are, of course, other precision constraints which one may take in to account. As 
we shall see in more detail, however, between B s — > Aa^, and b —> sj, the models 

we consider are already strongly constrained if we wish to use them to fit the experimental 



data within reasonable error estimates. We use the IsaTools 38-41 set of subroutines 



incorporated with ISAJET to calculate our estimates of these observables. 



IV. MINIMAL SUPERGRAVITY 



The first 
(mSUGRA) [42H46 



SUSY-breaking model we consider is minimal Supergravity mediation 
In the absence of other effects, gravity should act as a messenger 
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between the hidden sector where Supersymmetry is spontaneously broken and the Stan- 
dard Model sector. That is, the scale of the messenger interactions M mes is approximately 
the Planck mass Mpi, otherwise referred to as high-scale SUSY breaking. This leads to a 
gravitino/goldstino with mass on the order of a few TeV. Gravitational interactions induce 
SUSY- breaking terms at the high scale. 

In the minimal model, the GUT scale parameters are chosen to be a common scalar mass 
Mq, a common fermion mass M1/2, a common trilinear coupling A , and the value of tan/3. 
All other parameters are fixed except the sign of \i. We will consider only the fi > case, 
since otherwise we will have a larger than 3a discrepancy in Aa^ for the decoupling limit 
and worse for detectable superpartners. The GUT scale parameters are run down to the 
weak scale, resulting in a typical ratio of gaugino masses Mi : M 2 : M3 ~ 1 : 2 : 6. The 
lightest supersymmetric partner is typically a neutralino. 

In Fig. 1 we show 5a discovery contours (solid red) in the M ,Mi/ 2 plane for four choices 
of tan j3 in the mSUGRA model. We set A = and fi > 0. We present two contours, 
corresponding to 30 fb _1 of data running at 14 TeV and to 300 fb _1 at that energy. For 
the lower luminosity we apply the cuts as described above. For the high luminosity figure, 
we assume more restrictive triggers will be required to reduce the total recorded event rate 
to manageable levels. For tan/3 = 30 and higher, we raise the pr cuts to 150 GeV[47j and 
apply a reduced b-tagging efficiency = 0.5. We also modify the invariant mass selection 
so that only the leading two jets in px are considered as candidates for the Higgs mass peak. 
In our previous work we found that this strategy can improve the statistical significance and 
the signal to background ratio for high neutral higgs masses. For tan/3 = 20, these cuts 
are very restrictive and would offer little improvement over the 30 fb _1 results due to the 
relatively low masses accessible. We include a 300 fb _1 contour in the tan/3 = 20 frame 
using a pt > 75 GeV cut and e& = 0.5 with any pair of jets considered as a candidate for 
the Higgs decay. 

The dark gray regions are excluded for theoretical reasons such as tachyonic masses at 
the weak scale or lack of electroweak symmetry breaking. The solid blue region at low M\ i 2 
indicates charginos with mass M x + < 103 GeV, which have been excluded by experiment 
except in the case where M 2 > 1 TeV or sneutrino masses are less than ~ 200 GeV [34j. A 
more general bound including these cases can be put at M x + > 92 GeV. For high values of 
My 2 and relatively low M , the lightest slepton, a stau, becomes the LSP. This region is 
indicated by the solid light-gray area on the plot. 

The experimental value for Aa M is shown by the solid cyan line. The shaded cyan (for- 
ward slant hatched) region around it indicates a ±2cr error around it. We represent the 
experimental value and a 2a error for b — » 57 with a green dash-dot-dot line and yellow 
(backward slant hatched) shading. 2 Note that for this choice of model the measured value 
(solid yellow) does not appear; the edge of the shaded region shows the lower 2a limit. This 
is because SUSY contributions generate a negative correction to the SM value, while the 
experimentally measured value is slightly above the SM prediction. For B s — > fi + fi~ we have 
drawn dashed magenta contours for BF(B S — > fi + fi~) = 1. x 10 -8 and 5. x 10~ 9 . The higher 
value, corresponding to the smaller area on the plot, is approximately the current exclusion 
limit, set by CMS and LHCb. It should be noted that this is also roughly the value which 
provides the best fit for CDFs reported excess. The outermost line is not yet excluded by 
any experiment but should be reached by LHCb in the near future. LHCb should be able to 



2 We take the la error to be the sum in quadrature of the experimental and theoretical errors quoted above. 
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approach the SM limit with a few fb _1 of data, which in principle would push the excluded 
region out to arbitrarily high SUSY masses, depending on the errors. 

Current LHC data already strongly constrains some areas of parameter space. We include 
an exclusion bound based on LHC searches for SUSY particles in events with jets and missing 



energy 48j, this limit includes 1.1 fb _1 of data. It appears as as a blue, dash-dot line on the 
figure. However, the available bound was calculated based on a scenario with tan /3 = 10 
and so should not be taken as definitive here. We also include a current exclusion bound 
for the mass of Aq as a function of tan /3, based on LHC searches for the decay H — > t + t~ . 
(49| This bound is shown by the dotted black line on the plot. At low tan/3 it does not 
significantly extend the excluded region beyond chargino searches, but at higher values a 
significant region corresponding to lighter pseudoscalar masses is already ruled out. It was 
calculated in an MSSM framework using the M# Max scenario and so may not exactly 
reflect the bounds in the specific SUSY-breaking models we consider. 

Several comments are in order. First, we see that the discovery curve is roughly a quarter 
ellipse in Mq and Mi/ 2 . This approximately tracks the shape of a contour of constant Ma- 
An increase in Mq or Myz tends to increase Ma due to differential running of the Higgs 
parameters. Recall that in SUSY theories, the mass of A is given to first order by 

„ Ml - Ml - Ml cos 2/3 . 

Ml = — ^ ^ ? -. 15 

A 2 cos 2/3 v ; 

Ml and Mfj 2 are equal (set by Mq) at the GUT scale but differ at the weak scale due to 
different Yukawa and trilinear couplings that affect their running. M and My 2 contribute 
to the size of these differences and a higher initial scale M can lead to larger differences after 
running. As we increase tan /3 we can detect heavier Higgs through the increased Yukawa 
coupling to 6s and the discovery contour moves to higher Mq and M\i 2 - For tan/3 ~ 20, 
this range is up toM , M 1 / 2 ~ 250 — 300 GeV with 30 fb _1 of integrated luminosity. For 
tan/3 ~ 50, this extends M 1/2 out to ~ 750 GeV and M ~ 1600 GeV. With high luminosity 
running this range is extended to Mq, My 2 ~ 350 — 400 at low tan /3 and up to My 2 ~ 1000 
GeV or M ~ 2100 GeV at tan ~ 50. 

As expected, the regions excluded by or soon to be explored by B s — > fi + fi~ also grow 
rapidly with tan/3. For comparison one should bear in mind that the B s — » contours 
shown should be reached with only a few fb" 1 of data. The 30 fb _1 reach for the 36 signal 
significantly exceeds the current bounds from B s for all cases shown. However, at high tan j3 
the B s search will rapidly begin to outperform the 36 range at high Myz and relatively low 
Mq. Conversely, B s — > fi + fi~ performance is limited for high M and low M1/2, even at high 
tan/3. The leading order terms in tan/3 are proportional to M x ±/M^ which becomes small 
in this region (5lj . 

The 2a band around the measured value of Aa M prefers relatively light masses since we 
wish to generate a significant non-zero value. This range gradually moves to higher mass 
values as we increase tan /3. For tan /3 > 40 the preferred range is almost entirely covered by 
a 36 Higgs search, while for lower values the lower 2a region eludes us. The measured value 
is covered by the 36 search for all cases shown. 

A significant tension in these models is generated by the 6 — > S7 prediction. The mea- 
sured value is slightly in excess of the SM prediction, while mSUGRA generically predicts a 
negative contribution to the decay rate. Thus only the lower 2a bound appears on the plots 
and satisfying this constraint favors the decoupling limit with very high masses. The new 
physics contributions to 6 — > 57 include Higgs-quark loops, which give a positive contribu- 
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FIG. 1: Discovery contours (solid red) for the 3 6-quark Higgs signal in mSUGRA with 30 (LL) 
and 300 fb -1 (HL). We have chosen Aq = and jj, > 0. Dark gray regions are excluded by 
theory. Light gray indicates a stau LSP. The dark blue area is ruled out by current chargino search 
limits. The experimental values with 2a errors are shown for Aa^ (cyan, forward-slant hatched) 
and b — > S7 (yellow, backward-slant hatched; central value in green dash-dot-dot). B s — > fi + fi~ 
limits of 1. x 10~ 8 (current, lighter) and 5. x 10 -9 (darker) are indicated by dashed magenta lines. 
Current LHC exclusion limits are shown for cfp — > t + t~ (dotted black) and SUSY searches in jets 
plus missing Et (dash-dot blue). 10 



tion, and chargino-stop loops which unfortunately give us a larger negative contribution in 
this case. 

The result is that if we wish to satisfy the measured b — > 57 and Aa M values within 2a 
errors on both then we are pushed to a small region of parameter space in mSUGRA with 
A and sgn(n) as chosen. This overlap lies along the stau LSP region with M1/2 > 500 
GeV and low M . It is well known that the stau coannihilation region, where the lightest 
neutralino and the stau are nearly degenerate, provides one of the viable explanations for 
the observed dark matter relic density. Statistical fit analyses (52-54] over the parameter 
space, taking into account a large number of constraints including dark matter density and 
searches, favor a similar region as indicated on our plots. In general, for this model, the Aa M 
constraint favors lighter SUSY mass parameters, while LHC data pulls us towards higher 
values. Large values of tan /3 are favored because they partially ameliorate this tension. 

The preferred area is covered by the 3b search with 300 fb _1 for tan/3 > 40 but becomes 
more difficult to cover with lower tan /3. It should, however, be well-explored by B s decay 
in the near future for moderate to high tan /3. Inclusive direct searches for SUSY particles 



also have the potential to rule out or favor this region with accumulating LHC data |55l . 156 
With 100 fb"\ the LHC is expected to probe mSUGRA space up to M 1/2 ~ 1400 GeV at 
low M and M 1/2 ~ 700 GeV at high M . 



V. ANOMALY MEDIATION 

If tree-level soft SUSY-breaking terms arising from supergravity are suppressed, there 
remain loop-level contributions arising from the superconformal anomaly [57|, |58|. Such 
suppression can happen in extra-dimensional models where SUSY breaking does not occur 
on the brane which includes the SM sector. These anomalies generate mass terms which 
depend on the renormalization group beta functions and a mass scale set by the gravitino, 
M 3/2 . 

M M = ^M 3/2 (16) 
9i 

M J = l ( ^ + l^ 2 (17) 
A y = -^M 3/2 . (18) 

The resulting SUSY spectrum at low scales is significantly distinct from that found in 
mSUGRA. In particular, the gaugino masses have the ratio M% : M 2 : M 3 ~ 2.7 : 1 : 7.1 
where the gluino mass term has the opposite sign compared to the other two. This results 
in the lightest neutralino being primarily a wino, with the lightest chargino and neutralino 
nearly degenerate. This has important consequences for decay phenomenology, with a long 
lived chargino. The purely anomaly generated terms are renormalization group invariant. 
However, this leads to an important problem: it predicts tachyonic masses for the sleptons. 
To ameliorate this problem several solutions have been proposed which can generate posi- 
tive mass contributions from, e.g., bulk terms, gauge- mediated terms, new Yukawa terms, 



or higher order effects [57|, 1591461 



The minimal anomaly mediated model, mAMSB, is a simple phenomenological model 
which assumes a universal additional mass term M , which is added to the anomaly generated 
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scalar terms at the GUT scale 62|, 63]. The addition of this new non-anomaly mediated term 
breaks the RG invariance of the soft-SUSY masses. The deviation from the pure-anomaly 
case for scalars depends at first order on the Yukawa couplings. Hence, for first and second 
generation scalars the mass terms can remain close to their GUT-scale, diagonal relations 
and this represents a possible solution to the SUSY flavor problem. The complete set of 
GUT-scale parameters for the mAMSB model can be given by Mo, M3/2, tan/3, sgn(fi). 

In Figure 2 we show the 3b discovery curves in the M , M 3 / 2 plane for the four values of 
tan /3, as in the mSUGRA case. The color and pattern coding are the same as used before. 
The tachyonic region occurs where Mq is not large enough to offset the negative contributions 
to scalar mass coming from M 3 / 2 . This region grows as tan/3 is increased. As before, the 
discovery curves roughly parallel curves of constant Ma, which exhibit interesting behavior 
for large tan/3. The mass depends on the splitting between M m and M H2 at the weak scale 
and for low to moderate tan /3 the curve behaves similarly to the mSUGRA case. Increasing 
M or M3/2 tends to increase Ma due to the differential running and the starting point at 
the GUT scale. However, as can be seen in the plot for tan /3 = 40, the situation becomes 
more complicated as the down- type Yukawa couplings increase. 

This is because, for values of tan /3 in this range, we approach an approximate Yukawa 
unification where the couplings of the up and down-type quarks are very similar at the GUT 
scale. This does not occur with mSUGRA assumptions even for high tan (3 where the weak 
scale coupling become similar. Threshold corrections at the SUSY scale are generically large 
for the b-quark, whereas they are required to be relatively small for good Yukawa unification 
64 . Following Wells and Tobe, the leading finite threshold corrections for the b-quark go 
as 

_ y MSSM_ y SM /sinp ^ -gfrM-g tan P ^A/itan/3 
b ~ yf M /sin/3 ~ 12vr 2 M| 32vr 2 M f ' 1 ' 

These corrections should be roughly 5% or less for good unification, but for TeV squarks 
each term can be ~ 50% in general. For the mSUGRA plots shown above A t = at the 
GUT scale and becomes negative at the weak scale, thus both terms above are of the same 
sign and the correction is generically large enough to spoil any unification. In mAMSB on 
the other hand, A t is generated by the beta function as shown above and is positive at the 
GUT and weak scales. It's size is also large enough at the weak scale to be comparable to 
the gluino mass. The result is that it can significantly offset the corrections proportional to 
03, rather than add to them. This cancellation is not necessarily so precise as to guarantee 
exact Yukawa unification, but the difference between y t and at the GUT scale is ~ 10% 
rather than ~ 100% as in mSUGRA. 

The result of this is that the running of Mm and M# 2 can become very similar in mAMSB 
for high tan/3 and secondary effects become important. At low M 3 / 2 , increasing M 3 / 2 tends 
to make the Yukawa unification better, generating less difference between the Higgs mass 
terms during the running from the GUT to the weak scale and thus producing a lighter 
Aq. However, at the GUT scale the mass terms have a common positive contribution from 
Mo and a negative contribution which scales as M 3 / 2 . Due to inexact Yukawa unification 
and contributions to Mhi from the tau coupling, the beta functions for Mhi and M#2 
retain some difference and larger M 3 / 2 increases this initial difference. Hence, even with 
very similar running the initial splitting becomes large enough that M 3 / 2 begins to increase 
Ma- When M 3 / 2 is large and Yukawa couplings are similar, Ma is not particularly sensitive 
to M since it does not contribute to any initial splitting at the high scale. One can see that 
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FIG. 2: Discovery contours (solid red) for the 3 6-quark Higgs signal in mAMSB with 30 (LL) 
and 300 fb _1 (HL). We set /i > 0. Dark gray regions are excluded by theory. The dark blue area 
is ruled out by current chargino search limits. The experimental values with 2a errors are shown 
for Aa^ (cyan, forward-slant hatched) and 6—7-57 (yellow, backward-slant hatched; central value 
in green dash-dot-dot). B s — > ij> + h~ limits of 1. x 10~ 8 (current, lighter) and 5. x 10 -9 (darker) 
are indicated by dashed magenta lines. Current LHC exclusion limits are shown for c/>° — > t + t~ 
(dotted black). 
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at the right edge of the curve Ma actually increases as Mq is lowered toward the tachyonic 
bound. This is because the lepton contributions to the initial splitting and the differential 
running can become important. For high M 3 / 2 they can be negative and their magnitude 
increases as the offset M decreases. 

For tan j3 = 50, a detectable Ma can correspond to very high values of both M and M 3 / 2 
because of the effects discussed above. On the other hand, as can be seen in the figure, 
the tachyonic region becomes very large, while another disallowed region, one without EW 
symmetry breaking, grows for high Mq and lower M3/2. This leaves only a narrow band of 
theoretically allowed parameter space. 

As in mSUGRA, Aa M prefers relatively light superpartners which sit in a band at low 
Mo and M 3 / 2 . This band is only moderately dependent on tan /3. The regions probed by 
B s — > grow quickly with tan /3 and current exclusion limits begin to become important 

for tan > 30. The shape of the B s contour echoes that of the mass and discovery curves. 
Similar behavior is also seen for the b — > sj allowed region. The situation with respect 
to b — > 37 is, however, notably different than in mSUGRA models. For mAMSB, the 
chargino-squark loops give a net positive contribution to the decay rate. The additional 
Higgs contributions remain positive as before so the total correction relative to the Standard 
Model is the right sign to account for the observed value. 

Qualitatively, this change in sign can be understood as a result of the A t term. If we 
consider b to s transitions with each quark coupling to one side of a chargino-squark loop, the 
dominant terms come from the Higgsino-like coupling for the bottom quark paired with the 
Wino or Higgsino coupling of the strange quark. Before mass diagonalization, a left handed 
quark couples to a left squark via the Wino or to a right squark via the Higgsino. Thus 
we have one term proportional to the Higgsino-Higgsino term in the chargino mass matrix 
times the left-right coupling in the squark mixing matrix. The other term is proportional to 
the Higgsino- Wino mixing times the left-left or right-right coupling of the squarks. Most of 
these remain qualitatively similar between the mSUGRA and mAMSB scenarios, however 
the off diagonal stop mixing term can change signs between the two. Recall that this term 
is generically given by 

Mlr cot /3. (20) 

For mSUGRA as shown above, we assumed Ao = at the GUT scale and renormalization 
effects drive it to negative values at the weak scale, so the mixing term is always negative 
for positive fi. In mAMSB, on the other hand, the A t term as given by Eq. (20) is positive 
at the GUT scale and remains positive and large enough at the weak scale to make Mm 
positive. The sign of the loop involving the Higgsino-Higgsino coupling is therefore changed 
while the other loops remain the same. As mentioned before, a generic feature of AMSB 
is a nearly pure Wino for the lightest chargino and, by extension, a pure Higgsino heavy 
chargino. Thus the loops which have changed sign tend to dominate over those which do 
not: the latter depend on the chargino mixing angle which is typically small. The end result 
is that the stop-chargino terms can give a moderate positive contribution to b — > 57, unlike 
mSUGRA. Examining the figure we see that b — > 37 now disfavors superpartners and Higgs 
which are too light because they would give too much positive correction to the Standard 
Model. For tan /3 ~ 20 — 30 we do somewhat better than mSUGRA. With relatively low 
Mq and high M 3 / 2 we can satisfy both Aa M and b — > 57 within one or two a error on either 
measurement. As we push to higher tan/3, the tension increases as b — > S7 requires higher 
masses while Aa M still favors lower values of the GUT scale parameters. Only two small 
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wedges remain at tan /3 = 40 where both predictions can satisfy the experimental numbers 
with 2a errors on each. These wedges can be excluded by B s — > fi + fi~ in the near future 
and are well within reach of the 3b search. 

Based on current LHC data, exclusion bounds for mAMSB have been extrapolated in 



Ref. 165J1. The authors find that mAMSB is excluded for M 3/2 < 40 TeV when M < 1 



TeV, for tan/3 = 10. 

At tan (3 ~ 20 even the 300 fb _1 search does not cover a significant region beyond that 
excluded by chargino mass bounds. However the detectable space grows quickly with tan/3. 
If tan/3 ~ 30 the 3b search probes a portion of the favored area with 300 fb _1 , while 
B s decay does not currently put any new constraints on the parameter space. As tan j3 
increases 30 fb _1 quickly becomes sufficient to explore the experimentally preferred areas. 
For tan ~ 50 our search can probe into the theoretically allowed band, but this solution is 
strongly disfavored by Aa M . Direct searches for gluino and squark cascade decays at LHC 
with 100 fb" 1 are expected to probe to M 3/2 ~ 140 TeV for low M and to M 3/2 ~ 100 TeV 
for high Mq 



VI. GAUGE MEDIATION 

A third proposal for SUSY-breaking involves additional matter which is charged under 



the SM gauge interactions [67H69| . If this matter also interacts with the hidden sector then 
it can act as the messenger for breaking. The minimal model, mGMSB, assumes N pairs of 
fundamental SU(5) GUT representations, 5 + 5. In general, the gravitino mass is given by 

M * = Tib (21) 

where F is the characteristic scale of SUSY breaking in the hidden sector. Since the super- 
partner masses are given schematically by M ~ F/M mes , the gravitino will be the LSP if 
the messenger mass, M mes is much below the Planck scale. This is the default assumption 
for mGMSB. The mass of the messenger 5-plets is an input to mGMSB models, but the 
superpartner masses only depend logarithmically on it. The observable SUSY masses are 
parametrized by A = ,f s , where F s = F/Ca- Cn is introduced since the scale of SUSY 
breaking in the visible sector may be lower than that in the hidden sector. The complete set 
of high scale parameters is then A, M mes , N, Cq, tan/3 and sgn(fj). The scalar and gaugino 
masses are given by 

Ml = 2^(5(1)^)' + c/^f + C/^Y) (22) 

M A = k a NA-^, (23) 

47T 

where k\ = 5/3, = ks = 1, and Y is the hypercharge. C 2 = 3/4 for weak doublets and zero 
for singlets, C 3 = 4/3 for squarks and zero otherwise. Changing the number of messenger 
fields N introduces a relative splitting between the gauginos and sparticles since the former 
scale as N and the latter as yN. Since the messenger scale can be well below the Planck 
scale and the presumed scale of flavor physics, and since the initial masses depend only on 
gauge couplings, GMSB is a potential solution to the SUSY flavor problem. 

In the figures below we set Cq — N — 1 and plot in A, M mes . The dark gray excluded 
region on the left comes from the requirement that M mes > A. As expected, varying M mes 
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has only a weak effect on most of the lines. For fixed A, increasing M mes has little effect on 
the overall scale of SUSY partners. This is seen in the contours for Aa M , which are nearly 
flat. The SUSY corrections to a M come from chargino-sneutrino loops and neutralino-smuon 
slepton loops. They scale as ~ tan 0M^/Mg USY and are thus insensitive to changing M„ 



71] . Increasing M mes does gradually tend to increase the pseudoscalar Higgs mass due to 
increased running effects. 

The curve for 6 — > 37 exhibits some interesting behavior which requires a bit of expla- 
nation. For mGMSB, the predicted corrections to 6 — > 57 are generically positive but not 
for the same reasons as in mAMSB. The A t terms are negative at the weak scale as in 
mSUGRA, so the contributions from chargino-squark loops are negative. However, they are 
comparatively small in mGMSB, so that the positive contributions from charged Higgs loops 
can dominate and the total correction is a small net positive. The chargino-stop terms are 
suppressed due to the super-GIM mechanism. That is, they depend on the squark mixing 
matrix which is unitary and in the limit of degenerate squark masses the various terms can- 
cel as required for unitarity. For low scale SUSY breaking, running effects are small and the 
degeneracy of scalars with the same quantum numbers at the messenger scale is only mildly 
broken. As seen in the figure, increasing M mes at low scales decreases the contribution to 
6—7-57 due to the increased Higgs masses and the spoiling of the squark degeneracy. No- 
tably, although the overall correction to the SM value can become negative for large M mes , 
it does not exceed the lower 2a bound on the measured b — > rate. 

For reasonably low M mes we can comfortably reconcile Aa M and b — > 57 over a wide range 
of tan 0. B s — > fi + fi~ limits do not currently put a strong constraint on these favored regions 
although at high tan/3 they should begin to cover this region with enough luminosity. The 
preferred region also indicates excellent prospects for a 36 detection of the heavy neutral 
Higgs if tan > 30. Much of the space would still be covered by this search even for 
somewhat lower tan/3. For a discussion of current mass limits see Ref. [72J. The exclusion 
depends on the mass of the NLSP, which for our choice of parameters may be a neutralino 
or, at high tan /3 a stau. This puts a bound on the gluino mass of approximately 650 GeV, 
which corresponds to A ~ 80 TeV. Direct searches for sparticle pair production at LHC with 
100 fb _1 in the mGMSB model are expected to probe to A ~ 600 TeV in the 77 + 
channel. [7C 



VII. MSUGRA REVISITED 

As discussed above, minimal supergravity appears to only marginally satisfy experimental 
constraints in a few limited regions of parameter space. This is due to the conflicting pull of 
the g — 2 excess, which favors lighter SUSY partners, and the measured 6—7-57 branching 
fraction, which receives the wrong sign contribution from chargino loops in mSUGRA and 
therefore favors heavier masses. However, we have seen that mAMSB models can change 
the sign of the chargino contribution and that this can be traced back to the significant 
positive sign of A t at the weak scale. This suggests that by choosing a large, positive A in 
mSUGRA we can induce a similar effect which may improve the prospects for this model. 

In the figure below we plot results for mSUGRA as before but with Aq = 1.5 TeV as our 
initial condition. This has several interesting effects on the weak scale observables. First 
we see that, as intended, the overall correction to 6 — > 57 can be made positive and the 
experimental central value is predicted (solid yellow line) for Mi/ 2 ~ 200 GeV with a weak 
dependence on M . The contours for discovery and for Aa^ are not dramatically effected. 
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FIG. 3: Discovery contours (solid red) for the 3 6-quark Higgs signal in mGMSB with 30 (LL) 
and 300 fb _1 (HL). We set N = Cq = 1 and \i > 0. Dark gray regions are excluded by theory. 
The dark blue area is ruled out by current chargino search limits. The experimental values with 
2cr errors are shown for Aa^ (cyan, forward-slant hatched) and 6—7-57 (yellow, backward-slant 
hatched; central value in green dash-dot-dot). B s — > limits of 1. x 10~ 8 (current, lighter) 

and 5. x 10~ 9 (darker) are indicated by dashed magenta lines. Current LHC exclusion limits are 
shown for <p° — > t + t~ (dotted black). 
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On the other hand, the tachyonic region for lower values of Mq and M1/2 is notably larger 
than for the Aq = case. The stau LSP region also grows somewhat. The result is that 
it becomes possible to satisfy the b — > S7 measured value while keeping the predicted Aa M 
within the 2a lower bound. 

The region for Aa M exceeding the measured value is largely excluded by the requirement 
that we avoid tachyons. However, for tan /3 = 50 a small area where Aa M falls exactly on the 
measured value remains. Moreover, the preferred value for b — > 57 runs through this patch, 
so it is possible to have both predictions very close to the measured numbers for Mq ~ 700 
GeV and Mi/ 2 — 200 GeV. One should note though, that the region within 2a error for 
b — > 57 extends over the entire plot. It also appears that the area of exact overlap for 6—7-57 
and Aa M as currently measured may already be excluded by the — > t + t~ searches. As 
M1/2 grows the chargino-squark corrections can again become negative but the SUSY effects 
diminish on the whole so we do not fall below the lower bound. 

Another significant effect of positive A is seen in the contours for B s — > fi + {i~. They are 
significantly reduced, such that only at tan /3 = 50 does a limit of BF < 5 x 10~ 9 begin to 
cover parameter space that isn't already excluded by theoretical requirements. Additionally, 
the shape of the contours is slightly more complicated than the single ellipsoid seen with 
Aq = 0. For tan (3 = 50 the contours include a region at low M1/2 which quickly fall offs for 
M1/2 ~ 300 GeV, then grows again to cover an ellipsoidal patch from M1/2 ~ 500 — 1400 GeV. 
This can be qualitatively understood as follows: The ISAJET calculation of the branching 
fraction is based on an effective Lagrangian for flavor changing couplings between strange 
and bottom quarks (4lj : 

-Jgf e// = ~D R f D QiH d + D R f D [a g M g + a u M u fj u + a w M w ]Q L H* u (24) 

where 

a g = -^-fiM-g, a u = -—^iiA t , a w = ^-^fiM 2 . (25) 

M g ,M u and M w are diagonal mass matrices that depend on loop functions arising from 
gluino-down squark, stop-higgsino, and wino-up squark graphs respectively. The effective 
FCNC interactions are proportional to a function xfc which depends on a sum over these 
three terms with appropriate mixing coefficients. (See reference for details.) For our regions 
of interest the wino-squark loops are not as important as the other contributions. At low 
M1/2 both the gluino-squark and higgsino-stop loops are of the same sign, corresponding to 
a positive A t at the weak scale and generating a significant enhancement of B s — > [i + fi~. 
As M1/2 increases, the sign of A t at the weak scale changes due to running effects. For 
M1/2 ~ 1000 GeV the higgsino-stop term dominates as A t becomes large and negative, 
leading to a detectable excess in the branching fraction until the relevant particles become 
too heavy. However, for intermediate values of M1/2 the A t terms are of the right order to 
cancel the other contributions, giving us the trough seen in the graph. 

In general, the discovery contours in this high A scenario are similar to the A case. One 
can see that they cover most of the experimentally preferred region for tan > 40. This 
scenario is quite interesting to us, since it shows the possibility, particularly at high tan (3, 
to satisfy both the b — > 57 and Aa M measurements in a region which is not particularly 
sensitive to B s — > fi + fi~ but should be well within the range of 3b (and related) searches. 



18 



tan p = 20 



tan p = 30 




FIG. 4: Discovery contours (solid red) for the 3 6-quark Higgs signal in mSUGRA with 30 (LL) 
and 300 fb" 1 (HL). We have chosen A = 1.5 TeV and /x > 0. Dark gray regions are excluded by 
theory. Light gray indicates a stau LSP. The dark blue area is ruled out by current chargino search 
limits. The experimental values with 2a errors are shown for Aa^ (cyan, forward-slant hatched) 
and b — > S7 (yellow, backward-slant hatched; central value in green dash-dot-dot). B s — > fi + fi^ 
limits of 1. x 10~ 8 (current, lighter) and 5. x 10 -9 (darker) are indicated by dashed magenta lines. 
Current LHC exclusion limits are shown for (fp — > t + t~ (dotted black). 
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VIII. CONCLUSIONS 



The search for supersymmetry and the details of electroweak symmetry breaking is cur- 
rently in an exciting phase with the LHC quickly accumulating data at 7 TeV and expected 
to increase energy in 2014. Updated constraints on precision measurements limit the param- 
eter space in many models and some will rapidly improve with LHC data. A key signature 
of supersymmetry is the extended Higgs sector with two additional neutral particles and a 
pair of charged Higgs, unlike the Standard Model. 

To limit the rather large number of parameters in the MSSM, we have considered three 
simple models of SUSY breaking: mSUGRA, mAMSB, and mGMSB. We have delineated 
the parameters space regions of these models where a 3 6-quark signature arising from bA°, 
bH° production should be visible at LHC with 30 or 300 fb _1 of integrated luminosity. Since 
the Higgs to 6-quark couplings become large at high tan/3 (a scenario favored by Yukawa 
coupling unification), prospects for such a Higgs search are most promising at large tan/3 
values. The 36 signal will be complementary to direct sparticle searches, and will provide 
information on the heavy Higgs sector of the model. 

The rare decay B s — > /i + fi~ is also highly sensitive to tan /3 and already strongly constrains 
new physics at very high tan values. Limits from this decay mode are expected to improve 
quickly, potentially excluding important regions of parameter space. At the same time, CDF 
has seen a weak signal near the current limit which would strongly suggest new physics if it 
were to be confirmed. 

Meanwhile, the measured 3a excess of a M suggests a positive sign for fi and relatively 
light SUSY particles. For mSUGRA, this is in tension with the experimental value of 
6 — > sj, which is somewhat above the SM prediction, while the theoretical prediction receives 
negative corrections from chargino-squark loops. This pulls one toward higher masses to 
minimize the corrections and only a small region of parameter space is left where both 
numbers can fall within 2a deviations. It is possible that a large, positive value for Aq at 
the GUT scale can ameliorate this tension. For moderate to high tan /3, the 3b search will 
probe into the allowed region, as will B s — > in the near future. 

Anomaly mediated models provide one way to improve the predictions, since they give 
a positive contribution to b — > 57 for a positive sign of fi. This solution favors moderate 
(< 40) values of tan/3 so as not to increase b — > 57 by too much and to avoid the rapidly 
growing constraints from B s — > The 3b search will cover most of the preferred space 

for scenarios with tan/3 > 30 but will require a very high integrated luminosity to probe 
lower values. 

Perhaps the most natural fit to experiment, among the models we consider, is mGMSB. 
In this case, relatively light messenger fields lead to degenerate squarks which results in a 
cancellation of the negative chargino-stop terms in b — > sj, leaving an overall small positive 
contribution. This allows one to fit both the Aa M and 6—7-57 measurements easily. If this 
is indeed a hint at the nature of SUSY-breaking then the 36 search is quite promising and 
A°/H° should be discoverable at the LHC for a large range of moderate to high tan/3. 
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